Charged Radial Infall for Spherical Central Bodies 

J. FranklirQ and F. Morton-Park 
Department of Physics, Reed College, Portland, Oregon 97202, USA 

Abstract 

A massive, charged, spherical central body can be neutralized by attracting particles of opposite 
charge. We calculate the time it takes to neutralize these central bodies using classical mechanics, 
special relativistic mechanics, and finally, the "forced" trajectories of general relativity. While we 
can compare the classical and (special) relativistic times, and find, predictably, that the special 
relativistic neutralization time is longer, a comparison of these times with the general relativistic 
result is not as directly possible. We offer the final calculation as a demonstration of dynamics in 
a general setting and in particular, the structural similarity of the (general) relativistic problem to 
the other cases. 
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I. INTRODUCTION 



The common lore, that one "does not expect any macroscopic body to possess a net 
charge" [I] is certainly compelling observationally. But it invites the question, especially 
relevant during black hole formation: "How long does it take for a charged spherical body 
to neutralize?" Here we answer that question in three different ways. We consider the 
problem of charge neutralization for "classical" Newton-Coulomb black holes, and compare 
the time required to neutralize a central spherical body to the time required for a (special) 
relativistic neutralization, and finally, a full general relativistic calculation. 

In all cases, we have a specific model in mind - take a central body of mass M carrying 
charge Q = Nq where q is some fundamental unit of charge (the charge of the constituent 
particles making up the central body, say) and N is a large number. The neutralization 
process consists of arranging N test charges, each carrying charge —q, in a sphere of radius 
R and allowing them to fall radially inward under the influence of both gravity and the 
Coulomb force. Note that we ignore the force on the charges associated with the shell itself 
(which, as we discuss at the end, changes the ratio of Coulomb-to-gravitational interaction 
by a factor of one half). We will take our model central body to be a black hole, so 
there is a natural radius tq associated with the end of the infall process, the Schwarzschild 
radius (in the un-charged case). This model ignores the interaction of the "test masses" 
forming the charged spherical shell, and since each test particle falls radially in, they all 
take the same time to reach the radius of the spherical body. So we have a minimum 
time-to-neutralization, with no angular momentum, and no inter-particle interaction. These 
simplifying assumptions mean that we can focus on a single test particle, and our work 
amounts to repeatedly finding the time it takes that particle to go from R (at rest) to ro in 
various settings of increasing physical sophistication. 



II. NEWTON-COULOMB BLACK HOLES 

Consider a sphere of radius r with mass M distributed uniformly throughout its volume. 
The sphere sets up a gravitational force on massive particles - we can define the escape 
speed for test particles as the minimum speed needed to start at the surface of the sphere 
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and end at spatial infinity (at rest). That is just the usual: 




GMm 



(1) 
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But this equation can, in a sense, be turned around - instead of the escape speed v, we can 
ask for what radius tq does light itself "just" escape the gravitational field of the central 
body? From ([I]), with v = c, we have: 



This classical result is well-known: John Michell, the British geologist discussed this special 
radius in 1783, and Laplace independently developed the same idea in 1796, it is given in the 
first edition of Exposition du Systeme du Monde (although he removed it from subsequent 
editions). More of the history of these early "dark stars" can be found in [TJ. 

The applicability of the argument that takes us from ([!]) to ^ relies on canceling the ms 
appearing on either side of (jl]) prior to solving for tq - otherwise, the massless photon returns 
= trivially. It is interesting, then, that the same fundamental length (associated with 
a given mass M) appears in general relativity, where r is the "event horizon" associated 
with the Schwarzschild spacetime. The physical interpretation is also the same - that is 
the radius at which not even light can escape the gravitational field of the central body. 
Both because of the classical argument, and its general relativistic relevance, we will take 
our classical central bodies to have radius r$ = 2 ~$ M , and that will define the impact radius 
for the test particles that we send in to bring about the neutralization. 

The problem we need to solve, then, is: Given a uniformly charged sphere of radius r 
carrying total charge Q = N q and mass M, how long does it take a test particle of mass m 
and charge — q to reach r$ if it starts from rest at Rl The setup is shown in Figure [TJ 



FIG. 1. A test mass m carrying charge — q is released from rest a distance R from a central sphere 
of mass M carrying total charge Q. How long does it take for the test mass to reach the surface of 
the central sphere, t-q? 
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A. Newtonian Result 



To answer this question in the classical setting, we could work directly from Newton's 
second law: 

GMm Qq 

mr = - (3) 

r z 4 7r eo 

and integrate twice to find r(t), then invert r(t*) = tq to find the time t* at which the test 
particle collides with the central body. It is slightly easier (and more generalizable for future 
calculations) to start from the Hamiltonian, representing the total energy of the system: 

1 . 2 GMm Qq 

E = -mr . 4 

2 r 4ne r 

Because we want to compare the charged case with the un-charged case, we define the force 
ratio (3 = Q % { Mm^ so that f 3 = 

recovers the Newton (with no Coulomb) result. Then (kh 

can be solved for r 2 : 

f 2 = — + 2GM (1 + 13) -. (5) 
m r 

Finally, we can render each term dimensionless by setting r = pro and t = a to for dimen- 
sionless p and a. Then 

,2 2£ 2 CM _ m 1 

P ,2 = , n2 + : n2 (1 + /3) ~ (6) 



where p' = Now r is the Schwarzschild radius, but to is some characteristic time that 

is up to us - motivated by the first term on the right in ([6| and the definition r$ = 2G J A , 
we set 7 s = c, and then: 

p' 2 = — 2 +(l + P) -■ (7) 
mc z p 

All that's left to do is encode the initial conditions: Let R = 777*0, then p(0) = rj, and 
p'(0) = 0. Conservation of energy demands that: 

= — 2 + (1 + 0) - — > ^ = - J m c 2 (1 + /3) -, (8) 
mc^ r/ 2 77 

so that <Q becomes 

p' 2 = (i + /3) [i-i]. (9) 

The goal is to solve ([9]) for p(cr), and then find a* such that p(cr*) = 1, corresponding to 
r(t*) = r . If we rewrite (|9| in terms of dp and (icr, 

<i<r = — dp (10) 
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we can integrate directly from a = to cr* on the left, p = rj to 1 on the right, giving an 
expression for a*: 



4 yjrj — 1 + 7] 7r + 2 7] tan 
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a* = V 1 '-==L ^ZVl. (ii) 

With this expression, we can compute the time it will take for a charged sphere to 
neutralize. However, the particular solution found here is not extensible to the general 
relativistic version of the neutralization problem, and we would like to use the same approach 
in both the classical and relativistic cases. From that point of view, we want a method that 
will allow us to determine a* absent an explicit expression (ugly or not). We use a Runge- 
Kutta routine with a root-finding procedure to solve ^ for a* such that p(a*) = 1 as a 
function of /3 |9J . The result for 77 = 20 (so we start at a distance R = 20 tq from the center 
of the central body) is shown in Figure |2} As is to be expected, the neutralization time 
decreases as the electrostatic-to-gravitational force ratio gets larger (we are increasing the 
total magnitude of the force on the test particles). For comparison and as a test of the 
method, the exact result (11) is plotted in Figure [2] as well. 
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FIG. 2. Time to neutralization for a classical Newton-Coulomb central body with a test particle 
starting at rj = 20. The black points are the numerically computed time, the blue curve is the plot 



of (11). 



As this is a classical calculation, the speed limit of special relativity will surely be violated 
at some point. In these units, the speed of light is 1, and it is clear from ^ that when: 
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the speed of the test particle will be greater than c. For = 0, we see that p' < 1 for all 
values of p, since p > 1 (when p = 1, we are at r = r , and the test particle has collided with 
the central body) and -^j < 1. That makes sense when you think about how the radius r 
was chosen, but if /3 > 0, it is possible that the inequality in (12) is achieved by values of 
p > 1, leading to particle speeds greater than that of light prior to collision with the central 
body. In order to address this issue, we re-run the analysis starting from the total energy 
associated with special relativity 



B. Special Relativity 



The equation replacing Q in the relativistic setting is: 

GMm Qq 



E 



m c 
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(13) 



where the rest energy of the test particle is now included in E (and the rest of relativistic 
dynamics is built in). In the dimensionless variables of the previous section, this equation 
reads: 
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and we can solve for p' 2 to construct the analogue of Q, 
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Once again imposing the initial condition: 
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so that 
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Notice that this equation enforces a speed limit of 1 as it should [10] • The neutralization 
times for the model problem are shown in red in Figure [3] - as the force ratio increases, 
the special relativistic version has neutralization times that are longer than the classical 
calculation. This makes sense, since classically, particles can travel faster than c, while the 



relativistic speed limit . . . limits particle speeds. When the attractive force is greater, a larger 
portion of the classical trajectory has particle speed greater than c, while the relativistic 
trajectory always has speed < c, leading to a longer time. 
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FIG. 3. In red, the relativistic time-to-neutralization for our model problem with ij = 20, and 
in blue, the original non-relativistic result from Figure [2} At some point, the relativistic speed 
saturates, a feature not shared by the classical trajectories - since the special relativistic particles 
must travel < c, they have an overall longer neutralization time. 

There is, of course, a natural limit in the case of special relativity - the test particles 
need to travel a distance R — r and the maximum speed at which they could do so would 
be c, leading to 



t* = (18) 



or, in our current units, a* — r\ — 1. 



The Newtonian potential used in (13) is a natural starting point for relativistic investiga- 
tion, but the special relativistic case presented here does not take into account the fact that 
Newtonian gravity, unlike E&M, is not consistent with special relativity. The point-source 
solution we have used solves Laplace's equation, and therefore implies that the gravita- 
tional force changes instantaneously with a change in the source properties. While not 
an immediate concern for this static configuration, it is a hallmark of bad behavior under 
Lorentz transformation (in E&M, the magnetic field saves the special relativistic character 
of Maxwell's equations, and there gous approximate force that can be introduced 

to perturbatively correct the Newtonian gravitational field [2]). Indeed, the reconciliation 
of gravity with special relativity leads to general relativity. 



III. RELATIVISTIC BLACK HOLES 



In general relativity (GR), Einstein's equations relate the field of the theory to sources 
just as Maxwell's equations do for E&M. There are deep structural differences between the 
two, but the story we tell ourselves: "Sources generate fields" is the same. The field of 
GR (the target functions, like the potentials V and A of E&M) is the metric of spacetime, 
and the sources are mass and moving mass (more generally, since mass and energy are 
equivalent, any form of energy or moving energy). Once Einstein's equations have been 
solved, the metric influences the motion of test particles, just as once you have solved for 
V and A, you can find the force on test particles from the Lorentz force law in E&M. The 
nature of that influence, for gravity, is different, but the end result is a Hamiltonian similar 



to (13), depending in a more complicated way on the location of the test particle. The 



total energy replacing Q and (13) in general relativity, for a massive central body carrying 
charge Q (the so-called Reissner-Nordstrom [5] solution to Einstein's equation) is: 
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(19) 
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As a check, you should try expanding the above in r and verify that you recover (13). Notice 



that we only introduced the Coulomb potential as "additional" energy - the first term in the 
Hamiltonian is associated with the rest energy, the kinetic energy, and a "new" relativistic 
energy (reducing to the Newtonian gravitational potential in the appropriate limit). 
Now using our usual dimensionless variables, we have: 

+ — , (3 
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The new parameter (3 is again a ratio of electrostatic to gravitational forcing, this time with 
the central body's parameters only. We know that Q = N q, and if we take M = Nm (i.e. 
the test particles are precisely the charged particles making up the central body), then we 
have (3 = (3, and we can proceed with our usual routine. The resulting time-to-neutralization, 
as a function of (3 is shown in green in Figure |4| 
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FIG. 4. In green, the time-to- neutralization for the Reissner- Nordstrom relativistic solution. In 
red and blue are our previous special relativistic and classical results for comparison. 

In order to generate a roughly comparable time-to-neutralization here, we must be careful 
to define the starting point. In general relativity, the distance between two points (length) 
is determined by the metric of spacetime, and this metric has different values based on our 
choice of /3. If we insist that the (radial) distance travelled by the test particles be rji—rjf = 19 
(rjf = 1 by definition of our problem, and we start the particles at rji = 20), then each value 
of (3 used to generate Figure [4] requires a different starting point. In particular, we want 
i = 19, for: 

where i is the proper length associated with zero temporal separation in the Reissner- 
Nordstrom spacetime (i.e. a length calculated instantaneously). We can use this equation 
to determine the label of our starting point, rj*, in the modified spacetime setting |llj . 



1-1 + 



P 4p 2 



dp, 



(21) 



IV. CONCLUSION 



We have presented a simple model for the neutralization of a charged central body and 
carried out both the classical and special relativistic analysis for a variety of charge-to-mass 
ratios. The time it takes to neutralize a central body using special relativistic dynamics 
is longer than the corresponding time calculated from Newton's second law since in the 
relativistic case, particles cannot travel faster than c. In addition, we set up and solved 
the analogous neutralization problem in general relativity, using the Reissner-Nordstrom 
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metric to generate equations of motion for a charged test particle in that setting. It is 
interesting that while the general relativity form of the problem involves a variety of new 
physical ideas, the actual execution of the calculation is remarkably similar to its classical 
and special relativistic counterparts. As long as our test particles start far away from the 
central body, the general relativistic limit will match the special relativistic one (the test 
particles stay away from the strong field region of the spacetime, where interpretation is not 
straightforward). In our case, starting 20 Schwarzschild radii away led to virtually identical 
results for special and general relativity. 

Our model for neutralization was the simplest available - more sophisticated approaches, 
like that taken in [3] where the change in charge of the central body as the neutralization 
occurs is considered, provide additional physical insight when the classical and relativistic 
cases are compared. In addition, we could consider the force associated with the neutralizing 
sphere itself - that changes what we mean by (3 in a particular model by introducing the 
repulsive shell force: \ A ^ 2 for each particle This exercise, on the classical side, would 
provide a good introduction to the problem of spherical collapse in general relativity, where 
that same "self-force" is interesting |5]. Indeed in the general relativistic case, collapse 
depends on the details of the dust elements falling into the central black hole - the self- 
gravity of these cannot be ignored, and is an integral part of the study of astrophysically 
relevant collapse. Even our modest toy problem becomes more involved in the full GR case. 

In addition to the above simplifications, where we ignore the interaction of the test parti- 
cles with each other, we have also left out the interaction of the test particles with themselves. 
There are electromagnetic radiative corrections that will tend to slow the particle's progress, 
as energy is lost to the radiation fields. These corrections can be calculated in the classical 
and relativistic cases using familiar electrodynamics notions (as in Chapter 11 of |6J). 

In our simplified setting, we can put some relatively concrete bounds on the neutralization 
time. Taking the special relativistic result, in the large /3 limit, we know that the minimum 
possible time is a* — r\ — 1, or, with units: 



where rj is the number of Schwarschild radii away we use as the starting point, and ro is 
the Schwarzschild radius of a neutral black hole [12]. The minimum time scales like r\ - we 
expect r to be a small number - if we take a supermassive black hole with a mass of, say, 




(22) 
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10 6 M o , then r ~ 10 6 km and the time it takes for light to travel 20 r is about a minute, 
instantaneous on the timescale of the universe. Of course, the extreme special relativistic 
limit supposed that (3 is large, we can imagine much longer times, up to the amount of time 
it takes for uncharged material to fall radially into the black hole (the /3 — limit). 

From our assumption that /3 — /3, in which the infalling matter is the same as that 
constituting the central body, we have restricted ourselves, in the general relativistic version 
of the problem. For the Reisnner-Nordstrom spacetime, we must have Q < M (in units 
where G — c — 1), so the maximum value of Q is M, giving (3 < 1, the neutralization is then 
dominated by the charge-free case. If we relax this requirement, then we can have (3 = 1, 
say (for a maximally charged central body), and we are still free to choose (3 = ^ iv q e Gm - 
taking, for example, the proton leads to j3 ~ 10 18 , and that will give near-maximal values 
for the special relativistic case. 

The spherical shell represents a pretty minimal neutralization model - we expect more 
than just falling straight in from our neutralizing particles. In all of these settings, a more 
realistic "cloud" of particles with varying initial velocities and positions would lead to a more 
realistic (and necessarily longer) neutralization time, and these can be computed numerically 
in all three of the cases studied here. The inclusion of a self-force and calculation of a 
consistent spacetime for astrophysical black holes requires full GR and is an involved task 
- our goal here, motivated by the astrophysical observation that macroscopic bodies are 
neutral, is to provide a simple model that can be studied relatively completely in the classical, 
special and general relativistic frameworks. 
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